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. We calculate the gravitational correction to the phase difference between neutrino 

I mass eigenstates. There is a number of papers, where this correction was calculated 
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for the spherically symmetric gravitational field described by the Schwarzschild met- 
ric. The results of these papers differ from each other. Our result is close to that of 
Bhattachary, Habib and Mottola (see ref. [3,6]) and differs only by a coefficient: our 



' correction is twice smaller. 
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1. INTRODUCTION 

An investigation of neutrino oscillations is of great interest for particle physics. Neutrino 
oscillations in vacuum are described well by the existing theory (see, for example, [1]). 
However, the oscillation pattern changes for the propagation in an external gravitational 
field. In the paper we consider neutrino oscillation in the gravitational field described by the 
Schwarzschild metric. 

The correction to the phase difference of neutrino mass eigenstates due to the gravitational 
field described by the Schwarzschild metric was calculated in papers [2-8]. The results 
obtained in these papers differ from each other. In this paper we make our own calculation. 

In the problem of oscillations the spin of a particle is not important, so further we consider 
neutrino as a scalar particle. We also assume that the energies of the mass eigenstates which 
interfere are equal. Arguments in favor of such an assumption can be found in papers [9-12]. 

In Section 2 we calculate in the eikonal approximation the phase difference of neutrino 
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mass eigenstates propagating in a radial direction in the gravitational field described by the 
Schwarzschild metric. In Section 3 we compare our result with results of other authors. 
We use the system of units where h — c—1. 



2. RADIAL PROPAGATION IN THE SCHWARZSCHILD METRIC 

We consider neutrino that was produced at the distance ta from the center of the star 
which creates the gravitational field and was detected at r^. Moreover, the source and the 
detector are on the line on which the center of the star is located. 

Gravitational field of a spherically symmetric massive body is described by the 
Schwarzschild metric (see [13], Eq. (100.14)). In this metric the square of the interval 
looks like: 

ds^ = g^.^dx^dx'' = (^1 - ^) dt^ ~ ~ ^) ' ~ 

where g^i, is a metric tensor, Vg — 2GM is the Schwarzschild radius (M is the mass of the 
star), dfl is the angular part of the metric which is not important for a radial propagation. 

The wave function of neutrino can be written in the form (subscript a corresponds to the 
a-th mass eigenstate): 

*„(x,i) = A„(x,i)e^*"(^'*). (2) 

Let us calculate neutrino phase $a(x, t) in the eikonal approximation. This approximation 

means that the variation of /c^^ = —d^a{'^,t)/dx^ and Aa(x, t) is small on the wave length. 
Following [13] (§ 53, § 101) and [14] (§ 8) we write the eikonal equation in the form 

where nia is the mass of a-th mass eigenstate. 

The equation (3) coincides exactly with the Hamilton-Jacobi equation for a classical 
particle which was solved for the Schwarzschild metric in [13] (§ 101). Thus we can use the 
result for $a(x, i) from [13]: 



^^^-Et + L^+ [ (eUI---^) --iml + ^){l-^) ) dr, (4) 
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where E is the energy defined at infinity, L is the angular momentum (we consider radial 
propagation so we put L = 0). 
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Performing Taylor expansion in rr? jE"^ up to the second term we get: 

rA 

rB 

— /rA,(i-i(i-^)-^(i-^)>. («) 
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Then for the phase difference of a-th and a'-th mass eigenstates we obtain 

^aa' = - $a = ^ {m ' r a) + {tb - T ^) - In ^, (7) 



where 5m\,^^ = m\ — rri^, and 5m^,^ = — m^,. 

The first term in the formula (7) is the well known expression for the phase difference for 
the propagation of neutrinos in vacuum. The second term does not depend on r^, so it is not 
due to a gravitational field. This term is the correction of the order of 5m^/ E^ to the first 
term in the flat space-time. Therefore the gravitational correction to the phase difference is 
described by the third term. 

There is a factor of (1 — r^/r)/8 in front of m^/ E^ in the formula (6). It explains why 
the second and the third terms in the formula (7) have the same coefficient 1/8 and opposite 
signs. Moreover, the coefficients 1/2 in front of the flrst term and 1/8 in front of the second 
and the third terms appear as coefficients of Taylor expansion of \/l — x in x. Therefore the 
coefficients in front of the second and the third terms are deflned unambiguously by the well 
known coefficient in front of the flrst term. 

Let us make some numerical estimates of the contribution of the gravitational correction to 
oscillations of solar neutrinos ^. For the Sun we take = 3 km, (r^ — r^) ^ tb — 1.5-10® km, 
E — 1 MeV. In order to get the upper limit of the gravitational correction we assume VA — Tg, 
Sm^,^ — 2m^5m^,„ fa 6.4 • 10"^ eV^ (since Sm^,^ — 5ml,^{ml + ml,) < 2m^5m^,„), where 
m = 2 eV is the present upper bound on the neutrino mass, 5m^,^ = 8 • 10~^ eV^ (see [16], 
p. 541). Then the phase difference is fs 3 • 10''' + 6 • 10~^ - 2 • 10~^^ (each number in 



^ In fact, the neutrino oscillation pattern changes because of the Mikheev-Smirnov-Wolfenstein effect (see 
[15]). But to estimate the order of the gravitational correction we do not take this effect into account. 
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this sum equals to the value of the corresponding term in the formula (7)). Since the third 
term is much smaller then the first and the second terms, the gravitational phase shift is not 
observable. 



Our result is close to that obtained in papers [3, 6] and differs only by a factor in the 
second and the third terms of the formula (7). The authors of these papers got 1/4 while we 
obtained 1/8. Most likely, there is a mistake in the Taylor expansion of the square root in 
the formula (14) in [3] (and in the formula (12) in [6]). In a conversation with us during the 
4-th Sakharov Conference on May 22, 2009 one of these authors, Emil Mottola, said that he 
admitted the possibility of such an error. 



In the papers [7, 8] the gravitational correction was calculated only in the order brv? jE. 
The following expression for the phase difference of neutrino mass eigenstates was obtained: 



The absence of the influence of gravity in the first order in j E is in agreement with our 
result (see the formula (7)). The next order in Sm^/E where gravity appears non-trivially 
was not considered in [7, 8]. 

But having obtained the absence of the influence of gravity in the first order in Sm^/E 
the authors say that gravity appears implicitly. To show this they rewrite global variables 
in terms of local variables. 

As local variables the authors of [7, 8] use local energy 
proper distance Lp between the points ta and tb'- 



3. COMPARISON WITH LITERATURE 



3.1. Papers [3] and [6] 



3.2. Papers [7] and [8] 




TB-rAl ■ 



(8) 




Here they keep only the first order in Vg in the second equality. 
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Then the formula (8) can be rewritten as follows 



' g 



2£;(i-)(rB) V 2Lp \ rA tb. 
Such a transformation from global to local variables can be done for any quantity. There- 
fore, the obtained apparent dependence on gravity is not specific for the physics of neutrino 
oscillations. 

3.3. Papers [2], [4] and [5] 

The result for the gravitational correction obtained in the papers [2, 4, 5] differs greatly 
from ours. The authors of these papers state that the nontrivial correction arises already in 
the terms of the order of Sm^/E. 
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